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EIT DATA PROCESSING SYSTEM AND METHOD 

This invention relates to the field of electrical impedance tomography (EIT), and in particular to 
Err data processing systems and methods. 

There are many industrial two-component flow applications where it is required to monitor the 
flow rate of one, or both, of the flowing components. In a two component flow where the 
continuous component is electrically conductive, and where there is a conductivity contrast 
between the continuous and discontinuous components, a technique known as dual-plane Electrical 
Impedance Tomography (EIT) has the potential to be used on-line to obtain accurate estimates of 
the local disperse phase volumetric flow rate, the mean disperse and continuous phase volume 
fractions and the distributions of the local axial, radial and angular velocity components of the 
disperse phase. 

The invention is concerned with online data acquisition and processing of measurements from two- 
phase flows, whose axial velocities can be up to lOms"^ with ±5% precision. To be able to measure 
two-phase flows at these velocities requires data acquisition and processing speeds of 1000 frames 
per second per dual-plane in order to genemte flow velocity maps and this has not been possible to 

date, using conventional techniques. 

•I 

The following prior art documents are incorporated herein by reference. 

DICKIN F, WANG M, Electrical resistance tomography for process applications, Meas. Sci. 
Technol. 7 (1996), 247-260 

WANG M, DICKIN FJ, BECK MS, Improved electrical impedance tomography data collection 
system and measurement protocols. Tomography Technique and Process Design and Operation 
ISBN1813122467 

Tn, (2000), TMS320C6202/B Fixed-Point Digital Signal Processors, Texas Instruments 
Incorporated 

WANG M., (1994), Electrical impedance tomography on conducting walled process vessels, PhD 
thesis, UMIST 

SUMMARY OF THE INVENTION 

It is an object of the present invention to provide systems and methods for multiphase flow 
measuring that are robust, non-invasive and can deliver on-line multi-dimensional flow 
information. 

According to a first aspect of the invention, there is provided an electrical impedance tomography 
(EIT) data processing system, for acquiring and processing data from two-phase flows, comprising 
a dual-plane sensor, a plurality of digital signal processing modules configured in a data pipeline 
processing arrangement and a plurality of data acquisition subsystems in communication with a 
first one of said digital signal processing modules. 
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Preferably, the data pipeline processing arrangement is capable of acquiring and/or processing 
1000 dual frames per second per dual-plane. 

5 Preferably, said dual-plane sensor comprises an electrode array in communication with one of said 
data acquisition subsystems. 

In a preferred embodiment, said data acquisition subsystems each include one or more of a voltage 
controlled current source, an equal-width pulse synthesiser, a synchronised digital demodulation 
1 0 unit and an over-zero switch. 

The voltage controlled current source preferably comprises a parallel structure of eight AD844 
chips or equivalents, configured as four pairs in which the two inverting inputs of each pair are 
cascaded together with a current-setting resistor. The four negative and four positive current 
15 outputs of the AD844 chips are summed together to provide total negative and positive current 
outputs. 

Preferably, the voltage controlled current source further comprises a DC restore facility in which a 
capacitor and a resistor connected to the non-inverting input of each AD844 chip restore DC 
20 components and cancel the DC offset at the current outputs. 

In a preferred embodiment, the voltage controlled current source includes a potentiometer for 
amplitude balancing. 

Preferably, said equal-width pulse synthesiser comprises a clock signal for triggering an address 
generator to continuously output addresses to a pre-programmed memory, the memory output 
being connected to a digital to analogue converter, the digital to analogue converter providing a 
staircase signal output, characterised in that different sampling rates are provided at different signal 
frequencies so that the staircase signal output has the same time step-lengdi at all frequencies. 

Preferably, the synchronised digital demodulation unit comprises 16 sets of programmable gain 
amplifiers each having an analogue to digital converter, a strobed First-In First-Out (FIFO) 
memory and control logic and wherein the signal output of the equal-width pulse synthesiser 
triggers said 16 analogue to digital qonverters to acquire measurement data in parallel at predefined 
intervals. 

Preferably, said over-zero switch comprises two multiplexers, each having 16 electrodes mounted 
in one sensing plane connected to their 16 output channels and two flip-flops controlled by a 
switching command fcom one of said DSP modules and, ideally, a D-type flip-flop for preventmg 
40 data being written prematurely to the second of said flip-flops. 

According to a second aspect of the invention, there is provided a method of acquiring and 
processing electrical impedance tomography data from two-phase flows using the system as 
claimed in any of the preceding claims. 
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According to a third aspect of the invention there is provided a recording medium having recorded 
- thereon computer implementable instructions for performing the method described in the preceding.^ 
paragraph. 

5 The system has been carefully designed- in a modular fashion and can consist of several data 
acquisition modules and computing modules. The data acquisition modules include a voltage 
controlled current source with a DC-restore circuit, an equal-width pulse synthesizer unit to 
produce sjnitfaetic waveforms for electric field excitation, synchronized switching, and sampling 
control, and a synchronized digital demodulation imit consisting of sixteen parallel measurement 

10 channels. The computing modules include a powerful digital signal processor (for example 
TMS320C6202B/6713) with an 1EEE1394 conununication interface. Several DSP modules can be 
pipelined for a series of tasks ranging firom measurement control to image reconstruction to flow 
velocity implementation. The DSP system can implement two image leeonstruction methods: the 
single-step SCG method and the back-projection method. An online-updating cross-correlation 

15 algorithm has been developed and is to be implemented on the system to calculate flow velocity. 
This technique is less time and memory consuming compared to direct cross-correlation 
techniques. The system architecture, detailed modular design, and data processing algorithms are 
introduced below. With a little optimization on hardware, a speed of data acquisition of 1 164 dual- 
jGrames (2.383 million data points) per second has been achieved with a root mean square value less 

20 than 0.8%@80kHz using a water phantom. 



. The very high speed of data acquisition and processing afforded by the system of the present 
invention is achieved as a result of one or more of the following: 



25 • reducing the time required for multiple-tasks, such as data collection, image reconstruction 

and flow velocity implementation by the use of the data pipeline process with multi DSP 
(see Figure 1); 

• the speed of signal coupling due to the transient time of the improved wide bandwidth 
30 coupling circuit with differential discharge chcuits, therefore significantly reducing the 

transient time (see the SDD, Figure 10); 

• the residual voltage in capacitance of coupling due to the electric charge being switched at 
the zero voltage point, or zero charged energy of the capacitance charged by the sinusoidal 

35 - current, significantly reducuig the transient time (shown m the OZS of Figure 1 1); 

• the data collection operation speed, improved by the parallel structure of the SDD of 
Figure 8; 

40 • the time required for transferring and storing data to memory, reduced by the logical 

controlled FIFO memory (see the SDD and Figures 3 and 8); 



• the limited sampling speed of the analogue to digital converter limits the application to 
higih frequency in conventional metiiods but this limitation was overcome in the present 
45 system by employing a rolling sampling method (see SDD and Figure 9d); 
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• the time limitation of the control operation of the micro-processor, improved by the use of 
the logical control (not DSP), see SDD and Figure 8; 

• improved data transfer fix>m data acquisition subsystems to DSP to empty FIFO memory 
5 for the next operation, realized by the use of a fast DSP-FIFO interface (Figure 1). 

Using the system and method of the present invention, properties of the flow can be measured 
accurately even when tiie flow is highly non-uniform. Dual-plane ETT has the added advantage of 
providing on-line voliune fraction and velocity distributions which can be used for the purposes of 
1 0 process monitoring and control. 

BRIEF DESCRIPTION OF TBE DRAWINGS 

Preferred embodiments of the present mvention will now be more particularly described, by way of 
15 example only, with reference to the accompanying drawings, in which: 

Figure 1 illustrates the system architecture; 

Figure 2 shows one possible system configuration of a high performance two-phase flow meter 
20 using the DSP module as a building block, (a) Single Read-Operation Timing, (b) Single Write- 
Operation Timing 

Figure 3 illustrates the architecture of Ihe EIT data acquisition subsystems (DAS); 
25 Figure 4 shows a simplified schematic (a) and an equivalent schematic (b) for AD844; 
Figure 5 shows a voltage controlled current source (VCCS); 
Figure 6 shows the structure of the Equal-Width Pulse Synthesizer (EWPS); 

30 

Figure 7 shows waveforms generated by the EWPS; 

Figure 8 shows the structure of the Synchronised Digital Demodulation (SDD); 

35 Figure 9 shows the configuration of the sampling modes, (a) 4 samplings, (b) 8 samplings, (c) 16 
samplings, (d) rolling-sample method; 

Figure 10 shows an AC-coupling network; 

40 Figure 1 1 shows the trigger circuit of the over zero switch (OZS); 

Figure 12 shows a typical curve of independent measurements obtained from a phantom filled with 
mains tap water; 

45 Figure 13(a) shows photographs of two balls dropping, one relatively heavy and one relatively 
light; 
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Figure 13(b) shows 2D EU images; 

Figure 13(c) shows 1000 stacked images along section A in Fig 13(b) where a relatively light ball 
is falling through; 

Figure 13(d). shows 1000 stacked images along section B in Fig 13(b) where a relatively heavy ball 
is &lling through; 

Figure 13(e) shows the variation of averaged conductivities at the pixels along the sections A and 
B, according to time elapsed; 

Figure 14 shows cross-correlation curves at six sampling pixels along the radius of sensing planes 
via data collection speeds; 

Figure 15 shows 2D flow velocity profiles from an oil-water axial flow (M4) with superficial 
velocity: Vwater=0.249m/s, Voil= 0,05m/s; and 

Figure 16 shows flow velocity profiles from an air-water swirling flow with superficial velocity: 
Vwater=0.39m/s, Vgas=0.067m/s. (a) 2D profile from the combination of angular and radial 
velocity; (b) 3D profiles from the combination of all 3 dimensional velocity profiles. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Throughout these documents, tiie following terms are used: 

"on-line measurement" — synchronised measuring and imaging on the process line. 

**two-phase flow" — a fluid comprising two components, wherein there are diflFerences in the 
electrical properties of the two components. 

"SCO method" - a sensitivity theorem based algorithm using a conjugate grading method, 
"back projection method" - a linear approximation algorithm. 

"data pipeline processing arrangement" ~ an arrangement in which each processor (DSP) in the 
pipeline performs a particular function on the data and then hands over to the next DSP in the 
pipeline. All of the processors operate on data simultaneously so that data can be processed much 
more quickly than would be possible with a conventional DSP arrangement. 

System Overview 

The measurement system, as shown in Figure 1, includes up to four Digital Signal Processor 
modules (DSP1-DSP4), two identical Data Acquisition Subsystems (DASl, DAS2) connected to 
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two electrode arrays (Electrode Array 1, Electrode Array 2), and a third Data Acquisitioii 
Subsystem (DAS3) collecting data from conductivity, pressure, temperature and other auxiliaiy 
sensors/transducers. All three DASs are controlled by system-contFoliing DSP module DSPl 
through a common interface OData Acquisition Interface). Each of DSP1-DSP4 can be configured 
S to conmiunicate one at a time with a remote PC. 

All of the DSP modules are preferably identical to reduce the design and debugging tasks. The 
TMS320C6202B/6713 DSP chip \Tl 2000] may be used for the DSP module. Any DSP module 
can commimicate with a remote PC via an IEEE1394 interface, shown in Figure 1, to make the 
10 functionality of the system flexible. The four DSP modules work on a data pipeline processing 
technique to increase the data-processing speed. 

They can perform different tasks on data in the pipeline, for example in the illustmted embodiment, 
DSPl controls data acquisition and processes raw data; DSP2 performs image reconstruction, 
15 DSP3 fuses the image data to obtain flow information by performing velocity calculations and 
DSP4 is provided to carry out other specified functions. 

Figure 2 is a block diagram showing how DSP modules are used to implement such a high 
performance system. DSPl acts as the system controller and controls the data acquisition process 
20 of the three data acquisition subsystems (DAS1-DAS3) which are only coimected to DSPl. DSPl 
then formats the data and transfers it to DSP2. The data are transmitted via McBSPO on DSPl and 
received via McBSPl on DSP2. Besides these tasks, DSPl is also responsible for the IBEE1394 
inter&ce communication. 

25 EIT data acquisition subsystems 

Two identical data acquisition subsystems (DASl, DAS2) are included in the system to collect data 
from two electrode arrays (a dual-plane sensor). A block diagram of a DAS module is shown inside 
the dashed-line box in Figure 3. The four modules of most importance in the DAS module are: a 
30 Voltage Controlled Current Source (VCCS), an Equal-Width Pulse Synthesizer (EWPS), a 
Synchronized Digital Demodulation (SDD) unit and an over-zero switch (OZS) in the control base. 

Voltage Controlled Current Source (VCCS) 

35 The AD844 from Analogue Devices has a wide bandwidth of 60 MHz and very high ou^ut slew 
rate of up to 2000 V/ms, It differs from a conventional op amp in that the signal inputs have 
radically different impedances. As schematics shown in Figure 4, its non-inverting input (Pin 3) 
presents the unusual high impedance while the inverting input (Pin 2) presents low impedance, Rm, 
of about SOfl, A current applied to ttie inverting input is transferred to a complementary pair of 

40 wideband unity-gain current mirrors which deliver the same current to an internal slewing (TZ) 
node (Pin 5). It can be seen that Pin 5 and the inputs behave as a good VCCS with high output 
impedance and wideband. The one drawback is its small transient current range of ±10mA in 
comparison with geneml ERT requirements. To overcome this problem, a parallel structure 
consisting of eight AD844s has been developed and shown in Figure 5. 

45 
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Two inverse-phased sin-wave voltage inputs are buffered and then connected to the non-inverting 
input of each AD844, where two capacitances C and two resistances r can perfectly restore the 
Direct-Current (DC) components generated in the circuit and then cancel fhe DC of&et at the 
current outputs. Eight AD844s constitute 4 pairs of two-output VCCSs. Two inverting inputs of 
5 each pair of AD844s are cascaded together with a current-setting resistor R. The positive and 
negative current outputs of each pan are summed together to form larger current outputs. The total 
output currents can be estimated using the following function: 

±IouT =±2xVsDj/(2xRiN +R) (Equation 1) 

10 

Compared to other VCCS designs, the advantage of the VCCS shown in Figure 5 is its excellent 
balance property. The ± I^^ have ideal 180** phase difference and ignorable DC biasing and tiieir 

amplitudes can be easily balanced with the potentiometer, Rp. The maximum output current 
amplitude of this cncuit is 30mA. Its output impedance is about 7S0k SI ||18pF. 

15 

Equal-Width Pulse Synthesizer (EWPS) 

The EWPS structure, shovm in Figure 6, is shnilar to that of the widely used Equal-Sample Signal 
Generator (ESSG). An auto-run clock signal triggers the address generator to continuously output 
20 addresses to a pre-programmed memory. The memory then sequentially outputs data to a 12-bit 
Digital-to-Analogue Converter (DAC). A low-pass filter then smoothes the staircase signal output 
from the DAC to the required sine-wave voltage signal. The frequency, phase and amplitude of the 
generated signal can be set by outside controlling unit - DSP. 

25 The difference between EWPS and ESSG is that the former has different sapiple numbers at 
different frequencies to make the staircase wave use the same time step-length at all frequencies, 
while the latter has the same number of samples at all frequencies so the staircase wave has 
different time step-lengths at dififerent frequencies. This diJBference has led to a significant reduction 
in the complexity of the filter design using EWPS as it only needs to set one comer-frequency for 

30 the filter for all signal frequencies while the ESSG needs to set different comer-frequencies for the 
filter for different signal frequencies. 

Considering the use of the same step-length, fhe samplmg pulse frequency /, is unchanged. 

/, =^0 or /, =1/, (Equation2) 

N 

35 Therefore, 

CEquation3) 

Hence, fhe lowest harmonics from the EWPS is N and the EWPS may use one low pass filter 
with a fixed comer frequency for all waveforms. The step-length used in tiie EWPS is 0.5 ^s which 
40 gives a sampling pulse frequency of 2.0MHz. The waveforms generated by EWPS are illustrated in 
Figure?. 
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The specially designed EWPS also provides a convenient way to produce synchronized signals and 
sample control signals that are necessary in controlling SDD. The synchronous signal (square wave) 
is used to control the start and tibie end of the sampling sequences and tiie sampling pulses are used 
to actually trigger ADC's conversions at specific times. The timing of both signals and the width of 
5 sampling pulses play a very important role in the achievement of SDD accuracy. 

Synchronized Digital Demodulator (SDD) 

The Synchronized Digital Demodulator (SDD) is a high speed, synchronized parallel digital 
10 demodulation system, which is composed of 16 sets of Programmable Gain Amplifiers (PGA) and 
12-bit ADC, two sets of Strobed Fh-st-Ih First-Out (FIFO) Memory and control logic. Figure 8 
shows the functional block diagram of the SDD. The sampling pulse signal fi-om EWPS triggers all 
sixteen ADCs to acquire measurement data in parallel at precisely defined intervals. Two sampling 
modes are available. 

15 

In Mode I, 4, 8, 16, or 32 samples can be selected in one sine, wave* cycle as the same as the 
conventional sampling method. ' 

Mode n uses a rolling-sample method. It acquires all samples over several sine-wave periods with 
20 one sample per period and a fixed phase shift between one sample pulse and the one after. 
Sampling Mode n is particularly designed for sampling at a high frequency due to the speed 
limitation of the ADC. 

The configuration of the two sampling modes is shown in Figure 9. The completion signal 
2S combined from ADC's conversion status signals will then trigger a control clock to sequentially 
push converted digital data into two FIFO memories. All operations are fully logical controlled 
without the involvement of the DSP, except the starting signal of each firame of measurement. 

The coupling circuits, as shown in Figure 10, are made of a number of voltage buffers and passive 
30 components to provide a wide bandwidth and &st settling ac-coupling interface, as well as a dc 
bypass to the dc current bias of the programmable gain amplifier (PGA). The input can be 
configured in either the single or the differential input mode. Only the differential mode is 
illustrated in Figure 10 for the use of the adjacent measurement strategy in electrical impedance 
tomography, where R1=R3 and Rl is much smaller than R2. Therefore, R2, CI dommate the 
35 bandwidth of the coupling, Rl and R3 provide a DC bypass to the DC current oflfeet of PGA but 
the responding speed is fast since the transient time is dominated by Rl or R3. 

Over Zero Switch (OZS) 

40 Conventional ERT systems use sinusoidal current sources to excite the sensing field. Since a 
controlling unit, such as a PC, DSP or microcomputer, directly controls the current switching 
between electrodes, the switching operation may occur randomly at any time during a sinusoidal 
cycle. When the current value at the time of the switching operation is non-zero, electrical charge 
may reside in the capacitances of the signal passing circuit, such as the double layer capacitance at 

45 the electrode-electrolyte interface, the coupling capacitance and the stray capacitance, where the 

8 



wo 2005/022138 



PCT/GB2004/003670 



double capacitance plays a major role in all capacitances. These may result in residual potentials 
presented between these cq>acitances and produce an additional contribution to next excitation and 
measiirement. These residual charges will require several sinusoidal cycles to die-out, which means 
the whole data acquisition system (DAS) has to wait for several cycles after each switching 
5 operation in order to perform next voltage measurement. 

The over-zero switch (OZS) aims to eliminate the charged residual voltage by controlling the 
switching operation of the current injection. In an ideal case, no charged potential resides in either 
the double layer capacitance or the coupling capacitance after switching off the current injection so 

10 the settling time is shorter with the data acquisition speed being increased. Since the charged 
voltage has a iU2 delay with respect to the current, the switch should operated at the time either a 
zero-voltage presents between the double-layer capacitance or a maximum-current through the 
capacitance. The current-maximum is actually applied to optunise the switching time or the 
relevant phase difference between the synchronised signal and the sinusoidal signal, as the charged 

1 5 voltage between the double4ayer capacitance is hardly measurable. 

Figure 11 is the circuit used in the improved system to realize the over-current-zero switching 
technique. Two inverse-phased currents 1+ & I- are connected to the common inputs of two 
multiplexers Ul & U2 respectively. 16 electrodes mounted in one sensing plane are connected to 

20 the 16 output channels (CH[0:15]) of each multiplexer. Channel selection is together controlled by 
DSP command and a synchronize signal, SYN^SIN. SYN__SINf is a square-wave signal, which has 
the same over-zero point with two current signals 1+ & When DSP send out a switching 
command, an 8-bits data (each multiplexer is controlled by 4-bits data), ED[0:7], is written into the 
first flip-flop, U3. This data is later written into the second flip-flop, U4, by the next rising edge of 

25 SYN_SIN, and then make the two multiplexers switch to other channels. A D-type flip-flop, U5 A, 
is used in this circuit to prevent writing data into U4 before the data has been set-up at the output of 
U3. 

30 . Test results 

A laboratory test vessel (a **phantom") fitted with 16 stainless steel electrodes and fiUed with mains 
tap water was used in one of the trials. In the experiment, all the PGAs in the SDD were equally set 
to a fixed gain value in order to simplify the test. The current excitation fi:equency was 80 kHz. 
35 Measurements were taken fi-om all electrodes during each current excitation/projection, which gave 
256 measurements. Four samples per period were taken. Therefore, each fi-ame data has a total of 
1024 measurements. A simplified digital matched-filter was applied to demodulate the data. A set 
of independent measurements were then extracted and averaged firom their mirror data based on the 
reciprocal theorem. 

40 

The top and middle curves in Figure 12 present two of measurements obtained from vessel filled 
with mains tap water. The bottom curve denotes the relative change between the two measurements 
to reflect the repeatability of these measurements. To demonstrate the high speed of the new EIT 
system, a test was carried out by dropping two non-conductive balls. The two balls have similar 
45 dimensions but different densities (radius, SG and conductivity of the two balls are as follows: 
Rl=1.74cm, R2=1.75cm, SG1=2.55, SG2=1.13 and ai,2=0). The two balls were dropped to the 
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bottom of the phantom (see Fig. 13a), arriving at different times. The time diJBFerence can be 
clearly visualized from the EIT images (Fig.lSb), 1000 stacked EU unages along the sections A 
and B (Fig.lSc and 13d) and the reconstructed conductivity changes at relevant pixels as shown in 
Fig.l3e. The relatively heavy ball dropped in section A arrived at the bottom earlier than the 
relatively light ball dropped at section B, It can even be seen from both the photographs and EIT 
images (Figs. 13d & e) that the heavy ball boimced back when it first knocked the bottom. 

The total number of data collection points depends on the sampling rate (the number of samples per 
signal period). With a fixed sampling rate of 320kHz, there are 1024 measurement points from one 
sensing plane with 16 electrodes for one frame image at the signal frequency 80kHz, or 2.048 
million data collection pomts from a dual-plane sensor for 2000 frames of images at the signal 
firequency 80kHz. After tiie optimization of tiie controlling unit of DSP C6202 modules at aspects 
of data collection procedure, memory management, mtemal data transferring and program coding, 
the data acquisition speed of the system increases significantiy. A data acquisition speed of 1 163.6 
dual-frames per second (d^s) at 80kHz (2327.2 frames per second in total) can be achieved. The 
data acquisition speeds at different signal firequencies are summarized below: 



Data acquisition speed (d,a.s.) at 
different signal frequencies (s.f.) 



s.f. (kHz) 


10 


20 


40 


80 


320 


d.a.s. (dips) 


150. 
0 


295. 
4 


590. 
8 


1163.6 


928. 
2 



Data Processing Algorithms 

Two image reconstruction algorithms have been coded for the DSP module. They are Sensitivity- 
map-based Back Projection (SBP) method and a simplified SCG method - Single-step SCG 
(SSCO) method. The SSCG method is based on Ifae SCG method but only iterated once from a 
homogenous conductivity distribution. The SBP method is the simplest and fastest. It seems that 
the SSCG has a better performance than the SBP method. 

Data format in floating point is adopted in most cases due to Ihe use of Ifae ill-conditioned 
sensitivity matrix for image reconstruction. However, a DSP module C6202 with fixed-point data 
format was used in an earlier version of the system which was relatively slow to process the data in 
floating-point format It only can reconstruct images with the SBP method at a speed of 189.0 
frames per second, which is much lower than the online data collection speed. Therefore, a 
floating-pomt DSP module, C671 3 is used in the system of the present invention to obtain a higher 
image reconstruction speed. This has achieved an image reconstruction rate of 577.4 frames per 
second. 

An online^updating cross-correlation algoritiun has been developed Assuming that the two EFT 
electrode arrays are installed parallel on the flow-pipe, the axial-flow velocity distribution can be 
estimated by a direct cross-correlation metixod, as given in Equation 4: 

f,fi(^^'-r>)f2('»)^ « = 0,1,...,(N-1) (Equation4) 
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where N is the sample length, and fi(m) & f2(m) are the m"" up-flow and down-flow images 
respectively. Eq.l can be simply implemented online by updating the Ri2(n) with the new k-th 
images, as described in the following: 

i2<«(/i) = /2<r^(/i)+/j(Jt-/i)/,(it), « = 0,1,- sCN-l) (Equations) 

where superscript (k) indicates that the corresponding value is of k-th sample. This implementation 
can greatly save calculation time and reduce memory size. The computation speed is subject to the 
number of pau^ of images that are estimated as a maximum sampling length. 

The minimum acquisition time, t, can be taken as twice the product of the minimum transit time of 
the fluid, T, and the fractional velocity discrimination, k. Rearrangement gives the following 

k = t/2 (Equation 6) 

Therefore, the speed of data collection plays a dominant role in flie use of cross-correlation 
algorithm. If the distance between two sensing planes is d = 0.1m then x = O.lm/lOms"^ = 0.01s. Iff 
= 1000 frames/s was achieved then t = 0.001s. Using Equation 6 gives a discriminatory precision of 
0.05 or 5%. The local velocity can be derived from the correlated frame number, n, at each imaging 
pixel, as given by Equation 7: 

c = fd/n (Equation 7) 

The importance of data collection speed is demonstrated in Figure 14. The legend numbers are the 
X- and y- coordinates for a vessel with a diameter 20. Due to the limited data collection speed 
(Figure 14(a)), the correlated frame numbers (in respect to the lowest cc-value) are small and 
thei^fore, they produce a poor resolution in flow velocity measurement A good resolution from the 
fast ^stem is clearly illustrated in Figure 14(b). 
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